Primordial black holes from the inflating curvaton 
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The primordial black hole (PBH) formation is studied in the light of the inflating curvaton. The 
typical scale of the PBH formation is determined by curvaton inflation, which may generate PBH 
with 10 14 g < Mpbh < 10 38 g when curvaton inflation gives the number of e-foldings 5 < N2 < 38. 
The non-Gaussianity of the inflating curvaton does not prevent the PBH formation. 
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I. INTRODUCTION 

The origin of the structure in the Universe is due to 
the primordial density perturbations that exist already 
when cosmological scales start entering the horizon after 
the final reheating. 

Inflation can generate the primary source of the cosmo- 
logical perturbation, which is usually the vacuum fluctu- 
ation of a light scalar field. In the original inflation model 
the primary perturbation of the light field (i.e, the infla- 
ton perturbation) is converted into the curvature pertur- 
bation at the same time when the perturbation leaves the 
horizon. 1 

However, later studies revealed that there are many 
other candidates for the mechanism of generating the 
curvature perturbations after inflation [HI G3 . The 
idea of these alternatives is very simple; given that there 
are many scalar fields in the particle physics model, there 
will be many fields displaced from their minima at the 
end of the primordial inflation epoch, which leads to the 
multi-component Universe after inflation where any kind 
of isocurvature perturbations can be created. Then the 
isocurvature perturbations can source the creation of the 
curvature perturbations after inflation. 

In this paper we consider inflating curvaton [91413} that 
can lead to the PBH formation. Generation of the cos- 
mic microwave background (CMB) perturbations is sepa- 
rated from the PBH formation, since unification of these 
perturbations requires highly model-dependent argument 
that is not suitable for the purpose of this paper. Pre- 
viously, such unification scenario has been considered for 
the running-mass model or the Type-Ill Hilltop inflation 



1 This mechanism has been used for the PBH formation in Ref. [ij, 
where the secondary inflation is considered for the creation of 
the perturbation. This must be distinguished from the curvaton 
mechanism. 



model jl4 H18[ . PBH can also be produced from passive 
fluctuations [19|. (See also Ref.[2(j|.) Our strategy is to 
separate PBH formation from the original inflation sce- 
nario, so that PBH formation is liberated from the CMB 
conditions. Note however that the seed perturbations are 
generated during primordial inflation. 

Although not mandatory, the secondary inflation may 
solve cosmological problems of some specific models. For 
instance: 

1. The primordial inflation model may require e- 
foldings less than 60 so that the model can ex- 
plain the observed CMB. The short inflation could 
be needed to explain the spectral index and its 
runnning, or the inflation could be short because 
of the fast (rapid) rolling. 

2. The model may predict unwanted relics after re- 
heating. They must be diluted. 

3. The model may not have elementary dark matter 
(DM) candidate. Then the dark matter of the Uni- 
verse must be explained by something other than 
the particles. 

The most optimistic scenario of the inflating curvaton is 
that all these "problems" are solved by the secondary 
curvaton inflation. 

In the standard curvaton mechanism, the crucial as- 
sumption is that there are two components in the Uni- 



verse, where one component scales like matter (p a 



oc a 



while the other is the radiation (p r 



The inflating 



curvaton uses slow-roll (or sometimes fast-roll |21|) field 
instead of the oscillating one, which is a natural gener- 
alization of the usual curvaton mechanism. We consider 
the slow or fast-rolling curvaton (the inflating curvaton 
that scales like p a cx a _3e ™) in the presence of the ra- 
diation. Then the curvature perturbation is expressed 
as [Id 
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Here e,„ is defined as 



Pa = -3e w Hpa-. (3) 

Unlike the usual curvaton mechanism, the constancy 
of C, a requires some debate in the practical calculation, 
since there is no doubt that C, a is time-dependent when 
the radiation density is comparable to the curvaton den- 
sity. However, after the beginning of the curvaton in- 
flation, the ratio of the radiation decreases rapidly, and 
e w will soon behave like constant. In that sense the ini- 
tial quantity of Ca must be calculated in the inflationary 
phase slightly after a few e-foldings, which however must 
not be later than the significant evolution of the curva- 
ture perturbation [{J, [12[. In this paper t — t- ln i is the 
time when the initial quantities are defined. Note that in 
the inflating curvaton model significant evolution occurs 
when e w <C 1 , while in the usual curvaton it occurs when 
e w = 1. The significant evolution is accompanied by the 
significant change of the parameter r CT , which will be de- 
fined later. The curvature perturbation cannot evolve 
after r a ~ 1. 

In this paper we are assuming that the effective cur- 
vaton potential is always quadratic (= ±^m 2 (7 2 ) so that 
the ratio 5a j a behaves like constant after the horizon 
exit. Thinking about the scale dependence of the pertur- 
bation, "the initial 5a j a measured at the horizon exit" 
depends on the corresponding scale if a is moving during 
the primordial inflation, and finally the distribution of 
8a I a becomes scale-dependent. For the curvaton mech- 
anism, small a at the horizon exit means large pertur- 
bation tx 8a I a on that scale. 2 If the effective curvaton 
mass is determined by the main component of the Uni- 
verse, one can expect small mass during the radiation 
dominated epoch (22[. With these simple assumptions 
we examine the PBH formation in the generalized curva- 
ton mechanism. 



II. SCALE-DEPENDENT PERTURBATION FOR 
THE INFLATING CURVATON 

To show an illustrative scenario we consider a new in- 
flation scenario or a Type-I Hilltop inflation [1J] for the 



2 When we were finalizing the manuscript we found a new pa- 
per |2.'I . in which the usual curvaton mechanism has been con- 
sidered for a similar PBH formation. The secondary inflation is 
the crucial difference, since it expands the PBH scale and usually 
makes PBH heavier than the usual curvaton scenario. 




During the primordial inflation 



Curvaton Inflation 



FIG. 1: The curvaton potential is flipped during the primor- 
dial inflation because of the Hubble-induced mass, while it 
is negative during the secondary (curvaton) inflation. The 
density of the curvaton is subdominant during the primordial 
inflation. 



inflating curvaton, in which the curvaton potential sat- 
isfies e CT = \Mp (y^r^j < 1 during the curvaton infla- 
tion 

Consider the effective potential for the inflating curva- 
ton; 



V(a) = V + \m a {H) V, 



(4) 



where m 2 < is possible. Here the effective mass is de- 
termined by the main component of the Universe. As 
is shown in FigJTJ we are assuming m 2 > during the 
primordial inflation while to 2 < during the curvaton in- 
flation. Hereafter we use the index "1" for the primordial 
inflation and "2" for the curvaton inflation, respectively. 

In this paper we are assuming quadratic potential for 
the curvaton so that a and 5a can share the identical 
equation of motion. At the same time the equation be- 
comes a linear differential equation and consequently the 
ratio 5a I a behaves like constant. We are avoiding any 
deviation from the quadratic assumption. 

Let us see why the ratio can bahave like constant when 
the quadratic potential is assumed for the curvaton. For 
the simplest example, consider the equations when m a 
and H are constants: 



d 2 a n ^ r da o,„, 
^ + ^- +m l(H)a 



0. 



(5) 



where the equation has the solution of the form a oc e at , 
and the coefficient a is given by 
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Defining (3 = a/H, we find 

a(t) ~ a(0)e-P Ht = a{0)e^ AN , (7) 

where AA" = Ht is the number of e-foldings spent af- 
ter the horizon exit. 3 For the curvaton perturbation 
(a = <j + 5a), the equation of motion after the hori- 
zon exit gives d3 ffofcO) - < 5cr /c (0)e-' 3AAr . We thus find 
5(Jk(t) / cr(t) ~ Sak(0)/a(0) because of the cancellation of 
the evolution factor. Although the mass flips its sign 
during the evolution, one can take the initial condition 
everytime at the beginning of the epoch to find similar 
evolution factor. Here the growing solution (/3 < 0) is 
possible when the mass term has the negative sign [24j . 
The equation of motion gives either fast or slow-rolling 
solution when \r) a \ < 3/4 (i.e, when \m a \ < where 

2 

the eta-parameter is defined by r) a = jjfz- As a conse- 
quence of the cancellation, the ratio 5a /a behaves like 
constant as far as the evolution is given by the separable 
function. (See also Sec lII A\ ) 

For the inflating curvaton model, the source of the cur- 
vature perturbation is 5a, which leaves horizon during 
the primordial inflation and evolves thereafter until the 
curvaton mechanism starts to work. At the horizon exit, 
the spectrum is given by 

*..-(§)'. (8) 

where Hi is the Hubble parameter during the primordial 
inflation. 

Since we have two inflation stages (primordial and cur- 
vaton), we also need to consider horizon exit during the 
curvaton inflation. Indeed, Ref.pJ uses the perturbation 
that leaves horizon during the second inflation and ex- 
plains the small-scale perturbations. 4 . We are consid- 
ering the opposite scenario, in which the perturbation 
that exits horizon during the primordial inflation seeds 
the curvaton mechanism. The scale dependence of the 
perturbation is illustrated in FigJ3] We are focusing on 
the specific case in which the amplitude of the curva- 
ton perturbation becomes strongly scale dependent due 
to the Hubble-induced mass during the primordial infla- 
tion. For that reason we are not avoiding the 77-problem 
of the curvaton potential. 

In our scenario the seed perturbation (the curva- 
ton perturbation) created during primordial inflation is 



3 Here we have omitted the subscripts, but the evolution must be 
calculated for both inflationary stages. 

4 Another scenario has been considered in Ref. [14- 18] . Given that 
the curvature perturbation at the scale is given by 

the significant scale-dependence may appear either from t v {k) 
(running potential) or from Vs v {k). Our paper complements the 
approach given in Ref. |l7(l for the running e v (k). 



Vs a Horizon 




FIG. 2: The field perturbation has significant scale depen- 
dence. The peak appears at fcp BH , which corresponds to the 
scale that touches the horizon at the beginning of the curvaton 
inflation. In this picture the amplitude of the perturbation 
that leaves horizon during the first inflation decreases dur- 
ing primordial inflation (71 <S 1), while it increases during 
curvaton inflation (75 > 1). The perturbation entering into 
the horizon during the intermediate epoch causes oscillation 
around its average and decreases its amplitude. In the third 
picture, the perturbation newly created during the curvaton 
inflation is not dominating the amplitude at fcpBH- 

highly scale-dependent and it is converted into the curva- 
ture perturbation during the curvaton inflation. Our sce- 
nario complements the running-mass model or th e Ty pe- 
III Hilltop inflation model considered in Ref. |14T - ll8l |. 

For later use, we define the evolution of the perturba- 
tions as 

^ 2 (*ini)^ xTxTlnt%, (10) 

where t lri \ denotes the time when the inflating curvaton 
starts to work. Here the /c-dependent functions 71, Tint 
and Ti denote the translation functions during primor- 
dial inflation (Ti ~ e _/3lAJVl < 1 for ft > 0), interme- 
diate evolution (7^ nt ~ 1 for (3 int ~ 0) and the second 
inflation (T 2 ~ e -^ AN2 > 1 for < 0). Here AJVi 
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and A7V2 are the number of e-foldings elapsed during 
each inflation stage, defined for the perturbation on that 
scale (k = £;pbh)- AiVi is denned after the perturbation 
leaves horizon, and AN 2 is defined before inflating cur- 
vaton starts to work at t — t ml . T2 is usually negligible 
for the inflating curvaton. 

Here the inhomogeneity (5a) entering into the hori- 
zon can never be kept frozen even if the potential is flat. 
Oscillation around the average (not around the true vac- 
uum) starts just after the horizon entry and the ampli- 
tude is decreasing in the expanding Universe. If the os- 
cillation is sinusoidal, the evolution of the amplitude can 
be approximated as oc a~ 3 / 2 , which causes significant re- 
duction of the amplitude, as is illustrated in the second 
picture of Fig{2j We thus find a peak on the scale that 
touches the horizon at the beginning of the curvaton in- 
flation. We are using this scale for the PBH formation. 

For the PBH formation, we consider the spectrum of 
the source perturbation (Serf a), focusing on the scale cor- 
responding to the peak. The scale of the PBH satisfies 
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(11) 



Here Hi is the Hubble parameter during inflation, and 
pm is the energy density at the reheating after each in- 
flation stage. The reheating after the first inflation could 
be avoided, but for simplicity we are assuming the re- 
heating. 

Since the perturbation (fc = /cpbh) does not enter the 
horizon before the curvaton inflation, there is a bound 
for AiVi; 5 
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(12) 



where the scales are defined by Mf = "&H\Mp and M 2 = 
3iJ|M p 2 . 

Although the evolution of 5a is highly non-trivial, our 
final result does not depend explicitly on the evolution 
after horizon exit, as far as the curvaton potential is 



quadratic. That is because the ratio 5a j a behaves like 
constant after the horizon exit and thus the quantity 
evaluated just at the horizon exit is what we need for 
the PBH formation. Therefore, we find the perturbation 
using 5a / a, which is evaluated at the horizon exit |9lll2j: 



c 



2R 5a 



6t w a 

r a A 5a 
3r] 2 a ' 



(13) 



where the equation of motion is approximated by AH a ~ 
—dV/da with the coefficient A = (5 for the fast-rolling 
and A = 3 for the slow-rolling. Here r\ 2 denotes the value 
of r\ a during the curvaton inflation. The fraction 



3e u ,p a 



36^ ^ + 4/5R 
36.,,, 



3c™ + 4e- 



-4N 2 



(14) 



is defined at the end of curvaton inflation and throughout 

this paper we expect r a ~ 1 for simplicity. In the above 

2 2 

calculation we have used the ratio defined by R = 2 " Vo 
(m a < is assumed during the curvaton inflation) and 
the obvious relations 



p a = -3e w H 2 p 

dV . 
P ° " da-*- 



(15) 



The PBH is formed when the perturbation with the 
significant density contrast [5pj p > 0.1) enters the hori- 
zon. In our case the PBH formation occurs when the 
perturbation of the scale ~ /cpbh becomes accessible 
within the horizon. The typical mass of the PBH is given 

by mm 



k 



PBH 



0.00974Mpc~ 
which gives 

Mpbh — 



9* 

9*eq 



io« g f ll.) 
8 V100/ 



M, 



PBH 



6.67 x 10 50 g 



&PBH 



V lMpc~ 



(16) 



(17) 



where 3* is the degrees of freedom in the radiation. Let 
us see more details of the scenario. 

First, we need to check that £ at the PBH scale is 
dominated by the curvaton. Since the creation of 5a is 
possible also in the secondary inflation epoch, we need to 
compare it with the one generated during the primordial 
inflation. 6 We thus need the following condition at the 



5 The number of e-foldings during inflation is usually given by 
N = In °^, e i„fe^ ~ In a (f"? . In the above calculation the factor 

a{t)H\t) a(t) 

1/6 appears because H varies. 



The spectrum of &a, which is (newly) created during curvaton 
inflation, must not dominate the perturbation on the PBH scale. 
Since otherwise the curvaton mechanism is irrelevant for the PBH 
formation. 
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beginning of the secondary inflation; 

El <- Ele-PiAN! 

which gives a rather trivial condition 



M2 
Mi 



"1 

< e 2 



AJVi 



(18) 



(19) 



Second, since the PBH formation requires signihcant 
density contrast, we need 



V, 



1/2 r a A Hi 



3i] 2 2na kpBH 



0.1, 



(20) 



where 0fc PBH denotes the value of a when the perturbation 
(k = fcpBij) leaves the horizon. Then we find 



3r]2 J 2ir 



lO^GeV 
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Mi 



V3r/2 y V 1Ql5GeV 



(21) 



Third, £ at the CMB scale must not be dominated by 
the curvaton. The consistency at the CMB scale requires 



i/ 2 ^ rvA gi < ^ 

3/72 27TCr fccMB 

Therefore, from Eq.([20| and Eq.([22]l. we find 



(22) 



°"fcpBH ^ -p 1 (N 1 -AN 1 ) 



< 10~ 4 . 



(23) 



This condition suggests that a moves significantly during 
the primordial inflation stage. We thus find 



N x - ATVi > 9.2 x 
Here the scale can be given by 

&CMB 



PBH 



JV1-AW1 



(24) 



(25) 



where fccMB < IMpc . Note that steep potential allows 
significant running of the perturbation that leads to small 
k PBH (heavy PBH). 7 

According to Carr et al in Ref (26|, only PBHs with 
10 21 g < Mpbh < 10 28 g and 10 35 g < M PBH < 10 36 g 
can be allowed to become dark matter. 8 By observ- 
ing //-distortion due to the dissipation of large density 



7 Note however that Eq. H22t gives only the upper bound for the 
perturbation. The bound gives the heaviest PBH. 

8 CMB observations might have already exluded the latter pa- 
rameter space. In addition, quite recently it was reported that 
successful star formation could exclude PBH dark matter with 
10 16 g < Mpbh < 10 26 g H. 



fluctuation at a small scale after the decoupling of the 
double-Compton scattering, we will be able to check the 
amplitude of small-scale fluctuation by PIXIE |28| . 

As an illuminating example, consider Mi ~ 10 15 GeV, 
M 2 ~ 10 10 GeV and T m /M t ~ 10" 3 . These are giving 
a typical set of SUSY-GUT model. Then we find Ni - 
ANi ~ 10 and N2 ~ 30 (maximum curvaton inflation) 
gives AfpsH — 10 38 g, or N2 ~ 5 (minimum curvaton 
inflation) gives M PB h ^ 10 14 g. For A: PBH ~ 10 5 Mpc _1 , 
we find Ni - ANi ~ 11.5 and /3i > 0.8, which allows 
Mpbh - 10 36 g. 

In the original argument of the inflating curvaton 
model, there was a prediction for the running spectral 
index (n') Q. (Note that this "running" is not for the 
PBH perturbation. We are arguing here about the pos- 
sible relation between the PBH mass and the running 
spectral index of the CMB.) A similar argument can be 
applied for the PBH model, which suggests that large 
N 2 (i.e, small AT X ) gives both large Mpbh and enhanced 
running of the CMB spectrum. 

To be more precise about the relation between the 
CMB spectrum and the PBH mass, let us assume that 
the primordial chaotic inflation is driven by the potential 
oc <fi a and that it generates CMB through the conven- 
tional curvaton mechanism. In that specific model we 
find the spectral index n(k) — 1 ~ — 2ei = —a/2Ni(k) ~ 
-0.037 ± 0.014 with the running ri ~ -a/2N^ ~ 
-(0.037±0.014)/iVi. For the quadratic potential (a = 2), 
one will find Ni ~ 30 from the spectral index, which can- 
not be realized without the secondary inflation. The run- 
ning of the spectral index is enhanced because of small 
N\, which could be distinguished from other predictions. 
Note that in that model the second inflation is manda- 
tory and the condition ATj + N 2 ~ 60 - M p / H x ) 
gives the PBH mass M PB h < 10 3 V Here Ea.lfl^)). which 
bounds lS.Ni from below, is given by the scale of the in- 
flating curvaton and thus it can be related to the theory 
beyond the standard model. For instance, the gravity 
mediated supersymmetry (SUSY) breaking may predict 
the scale of the inflating curvaton at ~ 10 10 GeV and 
the inflating curvaton can dilute unwanted cosmological 
relics of the model such as gravitino or light moduli. 



For the fast-roll inflation [21| , the number of e-foldings 
is usually expected to be smaller than 60 {N\ < 60). 
In that sense the secondary inflation is mandatory for 
the scenario. Also, inflating curvaton could have a very 
interesting application to the rapid roll inflation [2j| . Us- 
ing the inflating curvaton, one can increase the inflation 
scale of this model to an extent such that gravity waves 
could be generated without worrying about curvature 
perturbation. Anyway, the curvaton-like mechanism is 
inevitable if the primordial inflation may not generate 
required CMB spectrum. If the primordial inflation is 
short, the mechanism must be accompanied by the sec- 
ondary inflation stage, which can be used cither for CMB 
or for short-scale perturbations. In that sense the inflat- 
ing curvaton is quite important for these models that do 
not have enough e-foldings by themselves. 
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There are another types of observation by which we 
can check the large running. If Vq is larger than > 10~ 6 
at small scales, Ultra compact Mini Halos (UCMHs) can 
be produced, which may be checked by observing future 
lensing events or gamma-ray due to dark matter anni- 
hilation [13, HH . The observation may help distinguish 
the running-mass or the Type-Ill Hilltop inflation model 
from the inflating curvaton model, since in the running- 
mass inflation model the perturbation on the CMB scale 
is fixed by the usual CMB observation and therefore its 
tail could be observed in smaller scales, while there is no 
such bound in the inflating curvaton model. 



A. Model-dependence 

Basically, the evolution of the perturbation before the 
curvaton mechanism is highly mode- dependent. In order 
to avoid the model-dependence, we considered the ratio 
da /a, for which the model-dependent evolution cancels 
when the quadratic assumption is valid. Any deviation 
from the quadratic assumption can lead to a highly mode- 
dependent result, which does not meet the purpose of this 
paper. 

In our calculation the ratio 5a /a evolves like constant 
before the beginning of the curvaton mechanism. Our ob- 
servation is very simple; since a and 8a are sharing the 
same equation of motion (when the potential is qudratic 
and the perturbation is beyond the horizon), and the 
equation is a linear differential equation, their ratio be- 
haves like constant. The typical evolution is shown in 
Eq.0 for constant m a and H, which gives the identi- 
cal evolution factor ~ e ~ l3AN for both Sak and a. The 
evolution is thus cancelled in the ratio. 



B. non-Gaussianity 

There have been many papers suggesting that the 
local-type non-Gaussianity could prevent PBH forma- 
tion [33| when Jml < —0.5. This bound could be crucial 
for the oscillation curvaton, since in that model one may 
find /nl — — 1 when r a ~ 1. On the other hand, the 
exact non-Gaussianity in the inflating curvaton scenario 
has been calculated in Ref.[12], which shows 



}nl 



1 5e w 

-5e ro -} 
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10 

i i 



10 
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(26) 



quadratic assumption we have 

/ = dg 
da* 



-9 



g'a* 



= 0. 



(27) 



where the ratio is defined as R = -k\m} a a 2 /Vq. Note 
that R is negative for the hilltop potential. With the 



Here the non-Gaussianity is evaluated for the PBH scale 
perturbation of the inflating curvaton. 

Unlike the oscillating curvaton, r a ~ 1 in the inflating 
curvaton does not lead to f^L ~ — 1- For r a ~ 1 there 
are cancellations in Eq. (|26]l . which finally gives Jnl > 
for R < 0. For r a < 1 there is no cancellation but Eq. (|2"6")l 
gives f NL > 0. 

We thus find that the sign of the non-Gaussianity pa- 
rameter is plausibly positive in our model. The problem 
of the oscillating curvaton is avoided in the inflating cur- 
vaton. 



III. CONCLUSION AND DISCUSSION 

In this paper we have shown that the PBH formation 
with the typical mass range 10 14 g < Mpbh < 10 38 g is 
possible if the curvaton potential is flipped due to the 
Hubble-induced mass. 

We have been avoiding highly model-dependent argu- 
ment, but the result suggests that the inflating curvaton 
can generate the PBH that could be the dark matter, 
and/or the PBH with the interesting scale that could be 
observable in future observations [13, [H[ . 

The secondary inflation stage, which has been consid- 
ered in this paper, sometimes plays important role in 
specific cosmological models. Although it depends on 
the inflationary model, the spectral index can be tuned 
by changing N\, which usually leads to enhanced run- 
ning of the spectral index. In other cases one may expect 
fast (rapid) rolling inflaton for the primordial inflation, 
which requires additional expansion. In both cases the 
secondary inflation is mandatory for the scenario. 

The inflating curvaton can also dilute unwanted cos- 
mological relics that are produced after primordial in- 
flation [32]. In the most optimistic case, the inflating 
curvaton solves above problems and at the same time it 
explains the dark matter of the Universe. 

Neither artificial interaction nor fine-tuning of the po- 
tential has been assumed for the mechanism. 

Our result also suggests that the Universe that has 
many isocurvature components may naturally lead to sig- 
nificant short-length perturbations. Even if they will not 
source PBH formation, they might leave observable sig- 
nal in the small scale perturbations. 

In addition, at small scales the second order induced 
gravitational wave signal can be larger than the first or- 
der in these types of running models [HI]. In Ref . [l8| . 
Alabidi et. al. studied this effect and showed that 
one could discriminate models by observing gravitational 
wave in future project BBO/DECIGO. Those observa- 
tions at small scales may complement CMB observations 
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at large scales, or UCMHs and /i-distortion at small scales 
to reveal the physics related to the PBH formation. 
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